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ABSTRACT 


Magnetic holes in the interplanetary medium are explained as 
stationary, non-propagating, equilibrium structures in which there are 
field-aligned enhancements of the plasma density and/or temperature. 
Magnetic anti-holes are considered to be associated with depressions 
in the plasma pressure. In this model, the observed changes in the 
magnetic field intensity and direction are due to diamagnetic currents 
that are carried by ions which drift in a sheath as the result of 
gradients in the magnetic field and in the plasma pressure within the 
sheath. The thickness of the sheaths that we consider is approximately 
a few ion Larmor radii. An electric field is normal to the magnetic field 
in the sheath. Solutions of Vlasov's equation and Maxwell's equations 
are presented which account for several types of magnetic holes, including 
"null-sheets", that have been observed. 



I. INTRODUCTION 


ORIGINAL PAGE IS 
OF POOR QUALITY 


In an analysis of high resolution (12.5 samples/s) magnetic field 

measurements obtained by Explorer A3 In the period March 18 to April 9, 

1971, Turner et al. (1977) found that most low-field Intensity regions 

4 

(B £ 1Y) were isolated and very thin (** 10 km or w 100 R^» where is 
the proton Larmor radius). They called these regions magnetic holes. 
Observations of such structures with two spacecraft (Fitzenreiter and 
Burlaga, 1977) show that the extent of magnetic holes is more than w lO^km. 
Magnetic holes form a class of current sheets. 

Three types of magnetic holes were identified by Turner £t al . (1977) and 
analyzed by Fitzenreiter and Burlaga (1977) : 1) holes with no change in the 

magnetic field direction across the current sheet; 2) holes with a 
reversal in magnetic field direction, but no rotation of the field vector 
across the current sheet; and 3) holes in which there is a rotation of 
the magnetic field vector across the current sheet. Turner et al . also 
exhibited current sheets in which the magnetic field intensity increased 
with no change in direction ("anti-holes”). The purpose of this paper is to 
apply the theory presented previously by Lemaire and Burlaga (1976) to explain 
the three types of magnetic holes and to explain anti-holes. 

For the calculations, we make the following assumptions: 1) magnetic 

holes are stationary inhomogeneities convected with the solar wind speed. 

2) They are bounded by thin current sheaths and the radius of curvature of 
these sheaths is very large compared to the thickness. As a consequence, 
we consider the geometry of the sheath to be planar, with a normal in the 

A 

z direction. 3) The physical variables are functions only of z. A) The 

A 

magnetic field,' £ (z), is normal to the z direction as for tangential 


1 


discontinuities. 5) In the frame of reference moving with the magnetic 

A 

hole there is an electric field jS (z) along z l.e., E « E - 0: this 

implies that there is no net plasma flow across the boundary layer. 

6) The plasma is quasi-neutral and consists of electrons and protons. 

Multi-ionic plasmas have been considered but will not be discussed in this 

paper. 7) The plasma is collisionless. 8) The diamagnetic currents are 

carried only by the protons (the electron velocity distribution is assumed 

to be isotropic). 9) The plasma is in thermal equilibrium on both sides 

of the current sheath and the bulk velocities at z * ±cc are both equal to 

zero in a frame of reference comoving with the solar wind speed. 

The basic equations are Maxwell's and Vlasov's equations. 

We obtained solutions of these equations which describe magnetic 

holes using the method described by Lemairc and Burlaga (1976). In all 

of the examples discussed below, we assume that the ambient magnetic 

intensity and particle density on each side of the current sheath are 

3 

5Y (nT) and 5 protons or electrons/cm respectively. For simplicity we 

also assume that the temperatures of the electrons and ions are the same 

4 

on both sides of the layer at z = ±co : T £ = T = 7.5 x 10 K; calculations 
with T £ f and with different temperatures on bot; sides have also been 
made but are not discussed in this paper. This complicates the structure 
of the layer but does not alter the basic physical effects, which are 
our principal interest for the present article. The following sections 
present solutions to the problems that we have posed above, and they 
include heuristic discussions, based on particle orbit theory, aimed 
at identifying the basic physical processes in magnetic holes. 
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A 

II, MAGNETIC HOLES WITH NO CHANGE IN B: Am » 0 

Let us consider magnetic holes in which the direction of magnetic 

* A 

field remains along x: JB * B(z) x. In particular, consider the case 

where B * IV (nT) in the middle of the hole, together with the other 

boundary conditions discussed in Section I. A self-consistent solution 

of Vlasov's equation and Maxwell's equations, obtained by the method of 

Lemaire and Burlaga (1976) is shown in Figure 1. 

The current in Figure 1 reaches a maximum at the two points where 
jdBj is largest. This follows directly from 


dB 


ti J-VxB-y -r— 
*o ~ — J dz 


* dB 

* dl- 


Note that J flows normal to B and z in this current sheet. jjJ is 
~ dB 


x 

in the middle of the magnetic hole because there — — goes to zero 

dz 

Jy changes from negative values to positive values. 

In the first orbit approximation is actually the sum of two 
currents, the gradient drift current 


zero 

as 


\ dB 
D “ B 2 dZ 

2 

(where W x * L H m v^ * Nk (T^+ T g ^) which results from the gradient drift 
of protons in the inhomogeneous field, and the magnetization current 


J 


M 


-__d 
d z 



(There is no current due to E, since electrons and protons drift in the 
same direction and with the same speed due to' ji x JB . ) Setting x 
( J D + J M ) * v x JL 8 ive s 


+ Nk (T; + T ) = constant, (1) 

— tx ex 

2 u 
^o 
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which show ii that the pressure is constant across the current sheet. In 
agreement with our assumption of equilibrium. The increase in N coincident 
with the decrease in B shown in Figure 1 is a consequence of the constant 
pressure condition, (1), since for the case considered the kinetic 
temperature is nearly constant across the sheath. 

In the example illustrated in Figure 1, we assumed that the kinetic 
pressure enhancement was due to a larger thermal proton and electron 
density inside the magnetic hole. But a similar structure would be 
obtained when the excess of particle pressure is due to a localized beam 
of suprathermal electrons or ions spiraling in an interplanetary magnetic 
flux tube. High intensity electrons spikes or protons bursts events 
sporadically detected in the interplanetary medium can produce such an 
excess kinetic pressure in an otherwise uniform solar wind plasma. 
Correlative studies would however be necessary to verify if some of the 
observed magnetic holes can be associated with these suprathermal particle 
beams . 

The electric field E in Figure 1 preserves local quasi-neutrality 
in the plasma. Its value is largest where is largest. This can be 

understood heuristically as follows. The electron velocity distribution 
is assumed to be isotropic. This 

implies that the bulk velocity of the electrons is zero and that their 
pressure tensor is isotropic with p g = NkT^. The momentum equation for 
the electrons is therefore given by 

dkT 
e 

dz 


- N e 


, „ dN ^ 4T 
kT — + N 
e dz 


dPe 

dz 
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which shows that the charge separation electric field is indeed 
proportional to |d Qn N/utj when dT^/dz “• 0; this is the case for the 
solution in Figure 1. If the temperature changed and N were everywhere 
constant, the same argument would apply and the jJEj needed to maintain 
equilibrium is this case would be greatest where j-~J was greatest. In 
general both N(z) mid T(z) can vary and they need not be in phase; in 
this case E z (z) could be complicated, but the physical processes would be 
basically the same as in the simpler cases Just described. 

The first inflection point in N(z) in Figure 1 is indicated by a 
dashed line; note the small accumulation of positive charge at a distance 
approximately IL 0 (the proton Larmor radius on side 2) to the right of 
this point and the small accumulation of negative charge at « Rj 1 to the 
left of this point. Since > B 

excess charge is closer to the inflection point of N than the peak of 
positive excess charge. This accounts for the asymmetry in the E z (z) 
profile in Figure 1. 

A charge separation An/N is required to support' JE. Poisson's equation 

dE ,2 

gives je| An » or An ^ indicating that there will be an 

dz dz“ 

accumulation of negative charge in the region where is decreasing most 
rapidly (i.e. , where the dN/dz is increasing with z) and an accumulation 
of positive charge where E^ is increasing most rapidly (i.e., where dN/dz 
is decreasing with z) . Such electrostatic double layers generally appear 
at the surface of a magnetic hole, having an excess kinetic plasma pressure. 
The excess of positive charges at the edges of a dense, plasma region 
results from the larger gyroradius of the ions compared to the gyroradius 
of the electrons. Indeed, the positive charges are deflected at a larger 
distance from the interior of the density (or pressure) enhancement than 


2 » Rq < R ^2 a ^d the peak of negative 
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the lighter negative charges. On the contrary an excess of negative 
charges is expected to accumulate at the edge of a plasma depression. 
Indeed, because of their small mass-to-charge ratio, the electrons tend 
to be reflected before the ions. This builds up a negative surface 
charge on the low density side and a positive charge accumulation toward 
the high density side (see Figure 2a). 

Willis U971, 1975) has discuss jd the charge separation electric 
field in a Ferraro-type magnetopause boundary layer, assuming that there 
is nc plasma in the magnetosphere (i.e. , on the low density side). In 
this extreme case when the density is strictly equal to zero outside the 
plasma region, only the positive charge accumulation is present at the 
edge of the enhancement (see Figure 2b). This "single-layer" structure 
is also described by the kinetic theory of Lemaire and Burlaga (1976) 
when the plasma density is forced to become equal to zero at z = +CG. 
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III. HOLES WITH A REVERSAL IN FIELD DIRECTION WITHOUT ROTATION (Auu - 180°) 

A 

Let us again consider the case where 1 JB * B(z) x, but now we allow 

B(z) to vary monotonically from B^ > 0 to B^ < 0. Clearly, B(z) must go 

through zero at some point in the magnetic hole in this case. Note that 

"neutral sheets" or "null-sheets" ideally correspond to this special 

case of magnetic field reversal without rotation. As before, let us 

3 

take B. = 5Y, B, = -5Y, N = N = 5/cm and T. * T *■ constant ■ 

1 Z 6 p l E 

4o 

7.5 x 10 K on each side of the current sheet. 

Figure 3 shows one posstble solution of Vlasov's equation and 
Maxwell's equations for the case cons:? lered. One sees that there is a 
single current peak, which is maximum in the center of the diamagnetic 
current layer, where |b] - 0. The density rises to a maximum at the 
center of the hole, as required by the pressure balance condition for 
equilibrium. The charge separation electric field is maximum in each 
of the two regions where N changes most rapidly, as required to maintain 
approximate local quasi-neutrality of the plasma. Each electric field 
region is maintained by a concentration of positive charge on one side 
of the maximum of |Ej and a concentration of negative charge on the other 
side; note that the positive charge concentrations overlap. Although 
the appearance of the magnetic hole in this case is rather different from 
that in Section II, the basic physical processes are the same in both 
cases. 

Figure 4 shows another solution of the Vlasov-Maxwell equations, 

A 

again for B = B (z) x and for the same boundary conditions that were 
used in the preceding example. Although Figures 3 and 4 are solutions 
for the same boundary value problem, the solution in Figure 4 is rather 
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different from that in Figure 3. This illustrates the well-known point 
that boundary conditions do not specify a uniciue solution to the 
Vlasov-Maxwell equations. The principal difference in this case is that 
Figure 4 shows a double layer current distribution while Figure 3 shows a 
single peaked current distribution. The current distribution in 
Figure 4 resembles that in Figure 1 in being double-peaked, and it 
resembles that in Figure 3 in having a non-zero current where B is a 
minimum. 
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IV. MAGNETIC HOLES WITH A ROTATION OF B 

Consider the case In which £ rotates 180° from one side of a 
magnetic hole to the other, and consider the same boundary ci 'it ions 
that were used In our previous examples. A solution of the Vlasov- 
Maxwell equations obtained by our numerical program is shown in Figure 5. 
The magnetic field intensity profile in this case differs qualitatively 
from those considered above in that there is now a secondary maximum in 
B where in previous examples B was a minimum. The reason for the 
secondary maximum is basically geometrical, as can be seen by considering 
the rotation of .B in the current sheet (see the insert in Figure 5) . 

When 1 JB rotates in a plane, the current has a component parallel to the 
local magnetic field direction as well as a component perpendicular to 

A 

B: 




dB 

3 

dz 




where B^ ■ B sin U? and B^ ■ B cos ip. Therefore if the magnetic field 

rotates without changing intensity (i.e., dB/dz » 0) the diamagnetic 

curreht is field aligned and proportional to the rate of rotation 

dat/dz. In Figure 5, the peaks in £ oc< lr near the points where the 

dB 

direction, uu, changes most rapidly (and where happens to be zero, 

as can be seen in the insert in Figure 5) . J is also non-zero at the 

dB 

secondary maximum of B even though -z- * 0 there, because the direction 

dz 


of JB is changing (~ f 0) at that point. 

The density profile in Figure 5 is anticorrelated with B(z) as 
required by the constant pressure condition. The electric field is 
again extremal at the extrema of dN/dz for the same reasons given 
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earlier. The charge separation profile shows a 

peak beside each extremum of E , as required to 

z 

charge peaks are separated fr the peaks in E 

z 


positive and 

maintain E , 
z 

by » 1 Rj, . 


negative 
and the 
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V. ANTIHOLES 
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Figure 6 shows a solution of the Vlasov-Maxwell ’ s equations for a linear 
"anti- hole" for the sane boundary conditions that we have used in previous 
examples. By definition, B is enhanced in some region with dimensions 

A 

on the order of several ion Larmor radii or more. Since B » B z, |jj 
clB 

is maximum where is extremum. We have assumed T ■ constant, so N is 
anticorrelated with B as required by the pressure balance condition. 
is again extremum where is extremum, and each peak in E is maintained 
by a charge separation over a distance « few R^. Thus, this type of 
anti-hole is physically equivalent to the type of magnetic hole described 
in Figure 1. 
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VI. CONCLUSIONS 


Magnetic holes are considered here as kinetic scale plasma 
inhomogeneities convected with the solar wind speed. These localized 
density and/or temperature enhancements are assumed to be sheet— like 
structures which are stationary in the solar wind flow. The plasma in 
the magnetic hole is separated from the background solar wind by sheaths whose 
thicknesses are only a few proton gyroradii. Diamagnetic currents 
flow in these sheaths. In the kinetic model presented here to describe 
the structure of these sheaths, we have assumed that the radius of 
curvature of the surface is mu r h larger than the thickness of the current 
sheath itself. Neglecting Coulomb collisions and other irreversible 
processes we have found, by using the models of Lemaire and Burlaga (1976), 
a variety of magnetic field structures similar to those observed in solar 
wind magnetic holes and anti-holes. Even null sheets , in which the magnetic 
field reverses with *|bJ going strictly through zero, can be modeled by 
the theory presented above (e.g., case III). 

There are a number of other factors that might be considered in 
constructing more refined, models of magnetic holes. Coulomb collisions 
although infrequent would tend to broaden the holes. In some cases the 
ion drift speed might exceed the Alfven speed and cause local instabilities. 
Variations in both ion density and temperature should be considered. 

The effects of alpha 

particle drifts and curvature of the current sheet might also be 
significant in some cases. These factors are essentially complications, 
however. The aim of this paper has been to explain the basic mechanisms 
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responsible for magnetic holes. It will be appropriate to consider the 
other factors when high resolution plasma observations are available 
which give n(z), T(z), etc.. In the magnetic holes. 
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FIGURE CAPTIONS 


Figure 1 Model of an isothermal magnetic hole in which [b| decreases 
to 1^ without changing direction and then increases 
symmetrically to its initial value. 

Figure 2 An illustration of the charge accumulation which produces 
an electric field in a magnetic hole and which determines 
the thickness of a magnetic hole. Case a) shows a low 
density on one side and a high density on the other side 
of a sheath. Case b) shows zero density on one side of the 
sheath. 

Figure 3 Mode of an isothermal magnetic hole in which the direction 

A 

of- B reverses without a rotation of B. 


Figure 4 Another isothermal magnetic hole in which the direction of 

A 

reverses without a rotation of B. Although the boundary 
conditions are the same as for the hole in Figure 3, the 
structures of the two holes are very different, illustrating 
the non-uniqueness of solutions of Vlasov's equation. 

Figure 5 Model of an isothermal magnetic hole in which the direction 

* o 

of B changes by a rotation of B through 180 . In this case 
there is a component of current along JB due to the rotation 
as well as a component of current normal to JB due to the gradients 

Figure 6 Model of an isothermal anti-hole. Compare with the hole in 


Figure 1. 
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